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An interesting method of obtaining equation of state for nuclear matter, from a density dependent 
M3Y interaction, by minimizing the energy per nucleon is described. The density dependence 
parameters of the interaction are obtained by reproducing the saturation energy per nucleon and 
the saturation density of spin and isospin symmetric cold infinite nuclear matter. The nuclear matter 
equation of state thus obtained is then used to calculate the pressure, the energy density, the nuclear 
incompressibility and the velocity of sound in nuclear medium. The results obtained are in good 
agreement with experimental data and provide a unified description of radioactivity, scattering and 
nuclear matter. 



I. INTRODUCTION 



Nuclear matter is an idealized system of nucleons interacting without Coulomb forces which is translationally 
invariant with a fixed ratio of neutrons to protons. The goal of nuclear matter theory is to obtain an equation of 
. state (EOS) for nuclear matter starting from the underlying two-body nucleon-nucleon (NN) interaction. The nuclear 
^ ' EOS, which is the energy per nucleon E/A = e of nuclear matter as a function of nucleonic density p, can then be 
used to obtain the bulk properties of nuclear matter such as the nuclear incompressibility, the energy density, the 
pressure, the velocity of sound in nuclear medium etc. The EOS is also of fundamental importance in the theories of 
nucleus-nucleus collisions at energies where the nuclear incompressibility K comes into play as well as in the theories 
of neutron stars and supernova explosions [1] and a widely used experimental method is the determination of the 
I incompressibility from the observed giant monopole resonances (GMR) [2]. In the present work we obtain an EOS 
. for nuclear matter using the M3Y-Reid-Elliott effective interaction supplemented by a zero range pseudo-potential 
representing the single-nucleon exchange term along with the density dependence which takes care of the higher order 
exchange effects and the Pauli blocking effects. The density dependence parameters of the interaction are obtained by 
reproducing the saturation energy per nucleon and the saturation density of spin and isospin symmetric cold infinite 
nuclear matter. One of the density dependence parameter is also utilized to provide estimate for the nuclear mean free 
path. The EOS, thus obtained, then used to calculate the pressure, the energy density, the nuclear incompressibility 
and the velocity of sound in nuclear matter. 

The M3Y interaction, which was derived by fitting its matrix elements in an oscillator basis to those elements of 
the G-matrix obtained with the Reid-Elliott soft-core NN interaction, has a profound physical basis. The ranges 
of the M3Y forces were chosen to ensure a long-range tail of the one-pion exchange potential as well as a short 
range repulsive part simulating the exchange of heavier mesons [3] . The real part of the nuclear interaction potential 
obtained by folding in the density distribution functions of two interacting nuclei with the density dependent M3Y 
effective interaction supplemented by a zero-range pseudo-potential (DDM3Y) was shown to provide good descriptions 
for medium and high energy a and heavy ion elastic scatterings [4,5]. The zero-range pseudo-potential represented 
the single-nucleon exchange term while the density dependence accounted for the higher order exchange effects and 
the Pauli blocking effects. Since the density dependence of the effective projectile- nucleon interaction was found to be 
fairly independent of the projectile [6], as long as the projectile-nucleus interaction was amenable to a single- folding 
prescription, the density dependent effects on the nucleon-nucleon interaction were factorized into a target term times 
a projectile term and used successfully in case of nuclear scattering [7] and cluster radioactivity [8]. Present work 
uses the same effective interaction within single folding model description to obtain nuclear EOS and provide unified 
descriptions of nuclear scattering, cluster radioactivity and nuclear matter. 
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II. THE DENSITY DEPENDENT REALISTIC EFFECTIVE NUCLEON-NUCLEON INTERACTION 



The general expression for the density dependent effective NN interaction potential v(s) is written as [7] 

v(s,p,e)=t M3Y ( S ,e)g(p,e) (1) 
where the realistic effective M3Y interaction potential [3] supplemented by a zero range pseudopotential t M3Y is given 

by 

t M3Y {s,e) = 7999exp(-4s)/(4s) - 2134cxp(-2.5s)/(2.5s) + J m (e)5(s) (2) 

where the zero-range pseudo-potential Joo(e) representing the single- nucleon exchange term is given by [5] 

Joo(e) = -276(1 - ae)(MeV.fm 3 ) (3) 

and the density dependent part [9] has been taken to be of a general form 

g(p,e)=C(l-p(e)p n ) (4) 

which takes care of the higher order exchange effects and the Pauli blocking effects. This density dependence changes 
sign at high densities which is of crucial importance in fulfilling the saturation condition as well as giving different 
Kq values with different values of n for the nuclear EOS. The value of the parameter n = 2/3 was originally taken 
by Myers in the single folding calculation [9]. In fact n = 2/3 also has a physical meaning because then j3 can be 
interpreted as an 'in medium' effective nucleon-nucleon interaction cross-section Co while the density dependent term 
represents interaction probability. This value of f3 along with nucleonic density of infinite nuclear matter po can also 
provide the nuclear mean free path A = l/(p <r ). Moreover, it also worked well in the double folding calculations 
with the factorized density dependence for the high energy heavy ion scattering [7] and the cluster radioactivity [8]. 



III. THEORETICAL FORMALISM FOR NUCLEAR MATTER CALCULATIONS 

The two parameters of eqn.(4), C and /3, have been fitted to reproduce the saturation conditions. It is worthwhile to 
mention here that due to attractive character of the M3Y forces the saturation condition for cold nuclear matter is not 
fulfilled. However, the realistic description of nuclear matter properties can be obtained with this density dependent 
M3Y effective interaction. Therefore, the density dependence parameters have been obtained by reproducing the 
saturation energy per nucleon and the saturation nucleonic density of the spin and isospin symmetric cold infinite 
nuclear matter. 

The energy per nucleon e obtained using the effective nucleon-nucleon interaction v(s) for the spin and isospin 
symmetric cold infinite nuclear matter is given by 

e = [3h 2 k 2 F /10m] + g(p, e)pJ v /2 (5) 

where m is the nucleonic mass equal to 938.91897 MeV/c 2 7 k F given by 

k% = 1.5tt 2 p, (6) 

is the Fermi momentum, p is the nucleonic density while p n being the saturation density for the spin and isospin 
symmetric cold infinite nuclear matter and J v represents the volume integral of the M3Y interaction supplemented 
by the zero-range pseudopotential having the form 

J v {e) = J J J t M3Y ( Sl e)d 3 s = 7999(4^/4 3 ) - 21 34(4^/2.5 3 ) + J 00 (e) (7) 

The eqn.(5) can be rewritten with the help of eqn.(4) as 

e = [3ft 2 4/10m] + [pJ v C(l - pp n )/2] (8) 
and differentiated with respect to p to yield equation 

de/dp = [fr 2 k 2 F /5mp] + J v C[l - (n + l)/3p n ]/2 (9) 
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The equilibrium density of the nuclear matter is determined from the saturation condition de/dp — 0. Then eqn.(8) 
and eqn.(9) with the saturation condition can be solved simultaneously for fixed values of the saturation energy per 
nucleon eo and the saturation density po of the spin and isospin symmetric cold infinite nuclear matter, henceforth will 
be called the standard nuclear matter, to obtain the values of the density dependence parameters and C. Density 
dependence parameters and C, thus obtained, can be given by 

= [(1 - p)Po n ]/[(3n + 1) - (n + l)p], (10) 
p={10me }/{h 2 k 2 Fo }, (11) 
k Fa = [1.57r 2 po] 1/3 , (12) 

C = -[2h 2 k 2 Fo }/[5mJ v p Q (l - (n + 1)M)L (13) 

respectively. It is quite obvious that the density dependence parameter obtained by this method depends only on 
the saturation energy per nucleon eo, the saturation density po and the index n of the density dependent part but not 
on the parameters of the M3Y interaction while the other density dependence parameter C depends on the parameters 
of the M3Y interaction also through the volume integral J v . The energy per nucleon can be rewritten as 

e = [3H 2 k 2 F /10m] - (p/p Q )[h 2 k 2 Fo (l - /?p»)]/[5m(l - (n + l)0 Po 1 )] (14) 

The pressure P and the energy density e of nuclear matter can be given by 

P = p 2 de/dp = [ph 2 k 2 F /5m] + p 2 J v C[l - (n + l)/?p n ]/2, (15) 

s = p(e + mc 2 ) = p[(3h 2 k 2 F /10m) + pJ v C(l - 0p n )/2 + mc 2 ], (16) 
respectively, and thus the velocity of sound v s in standard nuclear matter is given by 

v s /c = y/dP/de = \J [2pde/dp - h 2 k 2 F /15m - J v Cn(n + 1) p n+1 / 2] / [e + mc 2 + pde/dp] (17) 

The incompressibility K of the spin and isospin symmetric cold infinite nuclear matter which is defined as 

K Q - k 2 F d 2 e/0k 2 F - 9p 2 d 2 e/dp 2 \ p=Pa (18) 

can be theoretically obtained using eqn.(18), eqn.(6) and eqn.(9) as 

K Q = [-{3h 2 k 2 F j5m) - 9J v Cn(n + l)0p™ +1 /2] (19) 

Since the product J V C appears in the above equation, a cursory glance reveals that the incompressibility Kq depends 
only upon the saturation energy per nucleon eo, the saturation density po and the index n of the density dependent 
part of the interaction but not on the parameters of the M3Y interaction. 

IV. CALCUATIONS OF ENERGY PER NUCLEON, PRESSURE, ENERGY DENSITY AND VELOCITY 

OF SOUND 

The calculations have been performed using the values of the saturation density p = 0.16/to~ 3 and the saturation 
energy per nucleon e = — 16MeV for the standard nuclear matter. For a fixed value of 0, the parameters a and C can 
have any possible simultaneous values as determined from nuclear matter. Using the usual value of a = 0.005 /MeV 
for the parameter of energy dependence of the zero range pseudo-potential, the values obtained for the density 
dependence parameters C and have been presented in Table- 1 for different values of the parameter n along with 
the corresponding values of the incompressibility K - Smaller n values predict softer EOS while higher values predict 
stiffcr EOS. The form of C(\ — 0p n ) with n = 2/3 for the density dependence which is identical to that used for 
explaining the scattering and the a, cluster radioactivity phenomena [7,8] also agrees well with recent theoretical [10] 
and experimental [11] results for the nuclear incompressibility. 
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TABLE I. Compression modulus at different parameter values of n 



n 


P 


a 


C 


K Q 




fm 2 


MeV' 1 




MeV 


1/6 


1.0672 


0.005 


5.02 


201.2 


1/3 


1.1934 
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3.05 


236.2 


1 /9 


1 ^779 


u.uuo 


9 


971 ^ 


2/3 


1.6257 


0.005 


2.06 


306.3 


1 


2.3758 


0.005 


1.73 


376.4 


2 


9.1667 


0.005 


1.40 


586.8 


3 


41.4358 


0.005 


1.29 


797.1 
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In Tablc-2 the theoretical estimates of the pressure P and velocity of sound v s of standard nuclear matter have been 
listed as a function of nucleonic density e and energy density e using the usual value of 0.005/MeV for the parameter of 
energy dependence, given in eqn.(3), of the zero range pseudo-potential and also the standard value of the parameter 
n = 2/3. As for any other non-relativistic EOS, present EOS also suffers from supcrluminosity at very high densities. 
According to present calculations the velocity of sound becomes imaginary for p < 0.1/m" 3 and exceeds the velocity 
of light c at p > 5.2p and the EOS obtained using vu + TNI [12] also resulted in sound velocity becoming imaginary 
at same nuclear density and superluminous at about the same nuclear density. But in contrast, the incompressibility 
Kq of infinite nuclear matter for the Vu + TNI was 240 MeV while that by the present theoretical estimate is about 
300 MeV which is in excellent agreement with the experimental value of Kq = 300 ± 25 MeV obtained from the giant 
monopole resonance [13] and with the the recent experimental determination of Kq based upon the production of 
hard photons in heavy ion collisions which led to the experimental estimate of K = 290 ± 50 MeV [11]. 
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In Fig.-l plots of the theoretical estimates of the velocity of sound and the energy per nucleon e of the standard 
nuclear matter have been shown as a function of the normalized nuclear density p/ po using n = 2/3. The continuous 
line represents the energy per nucleon in MeV and the dash-dotted line represents the velocity of sound in units 
of 10~ 2 c for the standard nuclear matter. Fig. -2 presents the plots of the pressure P and the energy density e of 
the standard nuclear matter as a function of nuclconic density p using n — 2/3. The continuous line represents the 
pressure in MeVfm~ 3 and the dash-dotted line represents energy density in the same units for the standard nuclear 
matter. 



V. RESULTS AND DISCUSSIONS 



The theoretical estimate Kq of the incompressibility of spin-isospin symmmctric infinite nuclear matter obtained 
from present approach using DDM3Y is about 300MeV. The theoretical estimate of K from the refractive a-nucleus 
scattering is about 240 MeV-270 MeV [14,15] and that by infinite nuclear matter model (INM) [10] claims a well 
defined and stable value of Kq — 288 ± 20 MeV and present theoretical estimate is in reasonably close agreement with 
the value obtained by INM which rules out any values lower than 200 MeV. Using the value of the density dependence 
parameter j3 — 1.6257 fm 2 corresponding to the standard value of the parameter n = 2/3 along with the nucleonic 
density of 0.16/m~ 3 , the value obtained for the nuclear mean free path A is about 4/m which is in excellent agreement 
[16] with other theoretical estimates. 

There is, however, some uncertaity about the saturation energy per nucleon e = —lQMeV for the standard nuclear 
matter. A five parameter least square fit to the Bethe-Weizsacker mass formula for theoretical binding energies 
B(A,Z) Th = a v A - a s A 2 ' 3 - a c Z(Z - I) /A 1 ' 3 - a sym (A - 2Z) 2 /A + S where 5 = a p A-^ 2 for even N-even Z, 
= -a-pA- 1 ' 2 for odd N-odd Z, = for odd A, has been achieved by a 2 = {l/N)Y\B{A,Z) Th - B(A,Z) Ex } 2 
minimization, where B{A 1 Z)e x are the experimental binding energies from the recent Audi-Wapstra mass table 
[17] and the summation runs over all the N=3179 nuclei listed in the table which yields the new parameter values 
as a v = 15.25 ± 0.02 MeV, a s = 16.25 ± 0.06 MeV, a c = 0.6877 ± 0.0014 MeV, a sym = 22.18 ± 0.05 MeV 
and a p — 10.07 ± 0.85 MeV. Identifying a v as saturation energy per nucleon for standard nuclear matter and 
hence using the saturation energy per nucleon e = — 15.25MeV, the saturation density p — 0.16/m~ 3 , n = 2/3 
and a = 0.005/MeV the values obtained for the energy dependence parameters and nuclear incompressibility are 
j3 = 1.6149/m 2 , C= 2.01 and K = 295. IMeV respectively. 



VI. SUMMARY AND CONCLUSIONS 



In summary, we conclude that the present EOS which is obtained using the M3Y effective NN interaction, which 
was derived by fitting its matrix elements in an oscillator basis to those elements of the G-matrix obtained with 
the Rcid-Elliot soft-core NN interaction, has a profound theoretical standing. The form C(l — /?p 2 / 3 ) of the density 
dependence is also physical and provides excellent results for nuclear matter, a and heavy ion scatterings as well as 
a and heavy cluster radioactivities. The value obtained for the nuclear mean free path is in excellent agreement [16] 
with other theoretical estimates. The present theoretical estimate of nuclear incompressibility is in reasonably close 
agreement with other theoretical estimates obtained by INM [10] model, using the Seyler-Blanchard interaction [18] 
or the relativistic Brueckner-Hartree-Fock theory [19]. This value is also in excellent agreement with the experimental 
estimates from GMR [13] as well as determination based upon the production of hard photons in heavy ion collisions 
[11]. Thus, a unified description for nuclear scattering, radioactivity and EOS for standard nuclear matter has been 
achieved using a realistic effective interaction along with a physical density dependence. 
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FIG. 1. The energy per nucleon e and the velocity of sound v s in nuclear matter as a function of p/po- The continuous line 
represents the energy per nucleon in MeV and the dash-dotted line represents the velocity of sound in units of 10 -2 c for the 
standard nuclear matter using n = 2/3. 



9 



Pressure 




0.2 0.4 0.6 0.8 1.0 

p (frrr 3 ) 



FIG. 2. The pressure P and energy density e of nuclear matter as a function of nuleonic desity p. The continuous line 
represents the pressure in MeVfm~ 3 and the dash-dotted line represents energy density in the same units for the standard 
nuclear matter using n = 2/3. 
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